Spin-spiral inhomogeneity as the origin of ferroelectricity in orthorhombic manganites 
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We argue that the homogeneous distribution of spins in the spin-spiral state of orthorhombic 
manganites is deformed by the relativist ic spin-orbit interaction. The thus induced spin-spiral in- 
homogeneity gives rise to the ferroelectric response of the purely electronic origin in even-periodic 
magnetic structures. The mechanism is generic and explains the appearance of ferroelectricity in 
the twofold periodic structure of HoMnOs as well as the fourfold periodic structure of TbMnOs. 
Nevertheless, odd-periodic magnetic structures preserve the inversion symmetry and thus are not 
ferroelectric. Our analysis is based on the low-energy model derived from the first-principles elec- 
tronic structure calculations and the Berry-phase formalism for the electronic polarization. 



The possibility of switching the electric polarization P 
by means of magnetic field and the spin magnetization by 
means of electric field has attracted enormous attention 
to multiferroic materials because of their potential appli- 
cability in new electronic devises as well as the fundamen- 
tal interest to the problem of coupling between magnetic 
and electronic degrees of freedom [1]. The orthorhombic 
rare-earth manganites i?Mn03 (the space group Pbnm) 
are regarded as the key materials for understanding de- 
tails of such coupling. After discovering the switching 
phenomena in TbMnOs [2|, an unprecedented number 
of investigations has been carried out in order to clar- 
ify magnetic, structural, and ferroelectric properties of 
various i^MnOs compounds (e.g., Refs. [3|45|). 

All of them are improper ferroelectrics, where the in- 
version symmetry is broken due to some complex mag- 
netic ordering [1]. Since the Mn-atoms in the Pbnm 
structure occupy the inversion centers, the latter can be 
destroyed only if the magnetic unit is larger than the 
crystallographic one. There are two types of compounds 
with different periodicity, which are typically regarded 
as representatives of two main theories of ferroelectric- 
ity in RM1C1O3. The first one is (nearly) fourfold periodic 
TbMnOs, where the ferroelectric activity is believed to be 
due to the spin-orbit interaction (SOI) related spiral spin 
alignment [6]. Another one is twofold periodic HoMnOs, 
where P is typically ascribed to the (independent on the 
SOI) magnetostriction effect in the noncentrosymmetric 
E-type antiferromagnetic (AFM) structure [7[. 

Nevertheless, many questions remain. (1) There is no 
unique theory of multiferroicity in i?Mn03. Do HoMnOs 
and TbMnOs really behave as completely different sys- 
tems, where the ferroelectricity is caused by different 
microscopic mechanisms? Is it possible to unify these 
two cases? (2) When the temperature decreases, many 
spin-spiral manganites, including TbMnOs, exhibit a 
lock-in transition into a fourfold periodic commensurate 
phase [3] . What is the origin of these commensurability? 
Is it different from the twofold periodicity in HoMnOs? 
(3) The experimental magnetic structure, derived for 
TbMnOs and some related compounds by assuming the 
"spin-spiral" model, rises many questions: even at low 



temperature, in order to fit the experimental data, one 
had to rely on the elliptical deformation of the spin spi- 
ral [4|. However, in the elliptical distribution, the mag- 
netic moments at certain Mn-sites are substantially re- 
duced (up to about 3fiB [4), which clearly contradicts 
to the Hund's rule physics. Thus, the spin-spiral model 
is probably incomplete. If so, what is the true magnetic 
ground state of i?MnOs and how is it related to the fer- 
roelectric activity of these systems? 

In this work we will rationalize some of these questions. 
We will argue that the ground state of i?Mn03 is not the 
spin spiral. Instead, we will introduce the concept of the 
inhomogeneous spin-spiral, where the inhomogeneity is 
driven by the relativist ic SOI and is actually responsible 
for the ferroelectric activity in i^MnOs. We will argue 
that this concept is applicable to all even-periodic sys- 
tems, including TbMnOs and HoMnOs. 

The magnetic properties of manganites can be linked 
to the electronic structure of the Mn3d-bands located 
near the Fermi level. Thus, these bands, after the trans- 
formation to the real space, can serve as the Wannier- 
basis for an effective low-energy model. For i^MnOs this 
basis includes three ti g and two e g orbitals per spin for 
each of the four Mn-sites in the unit cell. The model 
Hamiltonian is taken in the Hubbard form, where all the 
parameters, such as the crystal field, SOI, transfer in- 
tegrals, and the effective Coulomb interactions are cal- 
culated rigorously by starting from the local-density ap- 
proximation (LDA) and using the experimental crystal 
structure [8] . The rare-earth 4/ states are treated as the 
core, which does not contribute to the low-energy proper- 
ties. The details of the computational procedure can be 
found in the review article [9] . The results of such calcula- 
tions for the whole series of the i?MnOs compounds with- 
out SOI have been reported in Ref. [10j. Here, we only 
emphasize two points, which are important for the mag- 
netic inversion symmetry breaking: (1) the Jahn- Teller 
distortion (JTD) gives rise to the large (~1.5 eV) crystal- 
field splitting between e g levels, which manifests itself in 
the orbital ordering (Fig. [I]); (2) The on-site Coulomb re- 
pulsion U is not particularly strong (~2.2 eV) due to the 
very efficient screening by the 02p band [10], which is im- 



FIG. 1. (Color online) Left panel: schematic view on the 
orbital ordering and magnetic interactions in the ab-plane of 
i?Mn03. Typical values of magnetic interactions for TbMnOs 
are Ji= -3.7, J-f= -0.2, J 2 b = -1.2, and J 3 = -2.3 meV 0, 
[Tl] ]. Rigt panel: dependence of the total energy (per one 
Mn atom) on the uniform spin-spiral vector q = (q x , q y , 1) 
obtained in the HF calculations for TbMnOs without SOI. 



portant for the "right" balance between nearest-neighbor 
(NN) and some longer-range (LR) magnetic interactions, 
whose form is controlled by the JTD [10|, 111 ). 

Besides JTD, the NN interactions in the ab-plane de- 
pend on the buckling of the Mn-O-Mn bonds. In the least 
distorted LaMn0 3 , Ji is ferromagnetic (FM) [To|, [H|. 
However, as the hybridization between t^g and e g states 
increases in the more distorted compounds, J\ can be- 
come AFM [10[. Since the on-site Coulomb repulsion 
is not particularly strong, besides conventional superex- 
change, there are other interactions, which appear in the 
higher orders of l/U, connect more remote sites, and 
compete with J\ . Among them, the third- neighbor inter- 
action J3, operating via unoccupied e g states of interme- 
diate Mn-sites, is the strongest one. The second-neighbor 
interactions Jf and J 2 are weaker. Nevertheless, what 
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important is the anisotropy \J^\ ^> | Jf I? which predeter- 
mines the direction of propagation (b) of the spin-spiral 
and the AFM E-state [10]. The NN interactions between 
adjacent ab-planes are strongly AFM (about —8 meV). 
The numerical values of J's for TbMnOs are listed in 
Fig. dJ Similar behavior was found for HoMnOs [101 ]. 

Without SOI, the competition of isotropic interactions 
in the ab-plane gives rise to the incommensurate spin- 
spiral state (Fig. [1]), which can be obtained by applying 
the generalized Bloch theorem [13] and solving the model 
in the Hartree-Fock (HF) approximation [14[ . Due to the 
AFM interactions between the planes, the spin-spiral vec- 
tor q, which specifies the phase cp^ = q • (r+R) + a T 
of the direction of spin = (cos(/?J^, suk/?^, 0) at the 
Mn-site r of the unit cell R, can be searched in the form 
Q = (#£C5 Qy: 1) 5 in units of reciprocal lattice translations. 
The total-energy minimum corresponds to the homoge- 
neous (a T =0) propagation along the orthorhombic b-axis 
with q y =0.6S for TbMn0 3 (and ^=0.72 for HoMn0 3 ), 
which exceeds the experimental values q y =0.2S and 0.25, 
reported for the be Q and ab helix, respectively. This 
discrepancy will be resolved later by considering the rel- 



FIG. 2. (Color online) Nonrelativistic electronic polarization 
for the general periodic q y — \ structure of TbMnOs depending 
on the angle <j) between spin magnetic moments in two Mn- 
sublattices. Relative directions of spins in the ab plane are 
explained in the inset. The numerical value of |P| in the E- 
state (0= 0° and 180°) of TbMnOs is 2.3 /xC/cm 2 (and 2.7 
/iC/cm 2 in the E-state of HoMn0 3 ). 



ativistic SOI, but first we turn to the analysis of the 
electronic polarization, which is believed to be closely 
related to the spin-spiral alignment [6|. For these pur- 
poses we compute P from the HF eigenvectors |C n k) in 
the Wannier-basis by using the Berry-phase formalism 
on the discrete grid of k-points [15[ . Since the Wannier- 
basis (and the model itself) was constructed by starting 
from the LDA bandstructure, which respects the inver- 
sion symmetry, it does not contribute to P. The details 
can be found in Ref. [l6[. First, we confirmed that with- 
out SOI, the homogeneous spin-spiral state does not pro- 
duce the electronic polarization. A finite P can be indeed 
obtained by switching on the relativistic SOI and per- 
forming one iteration by starting from the self-consistent 
nonrelativistic HF potential. As a test example, let us 
consider the q y = \ spiral in TbMnOs, which was inten- 
sively discussed in the literature [17j- Then, for the be 
and ab helix structure, we obtain P||c ~ 3 /iC/m 2 and 
P||a ~ 470 /iC/m 2 , respectively. Thus, we confirm that 
our minimal model successfully reproduces results of the 
first-principles LDA+U calculations for the spin-spiral 
state [17]: (1) the inequality P||a ^> P||c, which holds 
for the ab and be helix structures; (2) the absolute values 
of P, which are not particularly large. 

Being encouraged by this good agreement, we turn to 
the central part of our work, where we will show that (1) 
the homogeneous spin-spiral state in i?Mn03 is magnet- 
ically unstable as it tends to deform to some new inho- 
mogeneous state under the SOI; (2) this inhomogeneity 
gives rise to the new electronic contribution to P. 

The simplest example, which illustrates how the ferro- 
electricity can be induced by the spin-spiral inhomogene- 
ity without SOI is the general q y = \ periodic structure 
(Fig. [2]). In this case one can realize the collinear (and in- 
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FIG. 3. (Color online) Spin patterns in the ab plane of 
TbMnOs obtained in the HF calculations with SOI for q y — l 
(a), \ (b), | (c), \ (d), \ (e), and \ (f). The numbers 
indicate the angles between Mn-spins in different Mn-O-Mn 
bonds along the orthorhombic b axis. 



homogeneous in the sense that it is characterized by two 
magnetically different Mn-O-Mn bonds) E- state and the 
homogeneous spin-spiral state, depending on the angle (j) 
between two Mn-sublattices in the ab-plane: </>=0 and 
180° give rise to two E-phase domains with opposite po- 
larization, while 0=90° corresponds to the homogeneous 
spin-spiral structure with zero net polarization. Thus, 
in order to obtain finite P, it is sufficient to "perturb" 
the homogeneous spin-spiral state in the direction of the 
inhomogeneous E-state. For the collinear E-state this 
mechanism was proposed by Sergienko et at. [7[. Then, 
Picozzi et at. considered the behavior of polarization for 
the more general noncollinear alignment and on the basis 
of first-principles calculations argued that there is a large 
electronic contribution to P [18], which is semiquantita- 
tively reproduced by our model calculations in Fig. O 

Now the question is how to stabilize this inhomoge- 
neous spin-spiral structure? One possibility is of course 
the nonrelativistic exchange striction, which leads to the 
off-centrosymmetric atomic displacements and stabilizes 
the AFM E-state @,[l8|. We will come back to the anal- 
ysis of relative roles played by different mechanisms at 
the end of the paper, but first we want to show that the 
inhomogeneous magnetic structure, which gives rise to 
the finite ferroelectric polarization, can be naturally sta- 
bilized by the relativistic SOI. This mechanism is very 
generic and takes place in the twofold periodic HoMnOs 
as well as fourfold periodic TbMnOs. Since the general- 
ized Bloch theorem is no longer valid in the relativistic 
case, we consider the super cell geometries correspond- 
ing to q y =l/L for which we start from the homogeneous 
spin-spiral state, turn on SOI and further iterate the 
HF equations until self-consistency. Typically, this pro- 
cedure requires several tens of thousands of iterations, 
which are accompanied by the decrease of the total en- 
ergy. The results of such calculations for TbMnOs are 
summarized in Fig. [3l The equilibrium magnetic struc- 



FIG. 4. Left panel: electronic polarization in TbMnOs de- 
pending on the size of the magnetic unit cell L. Right panel: 
dependence of the total energy on L—l/q y for the homoge- 
neous spin-spiral states of TbMnOs without SOI (filled cir- 
cles), stabilization energy (AE) caused by the SOI-induced 
magnetic relaxation (open circles) , and total energies for mag- 
netic superstructures shown in Fig. [3] (squares). 



tures are characterized by an alternation of the angles 
formed by spins in different Mn-O-Mn bonds, which indi- 
cates the inhomogeneity. The obtained spin patterns re- 
flect the competition of many interactions in the system: 
the isotropic interactions tend to form the uniform spin 
spiral, while relativistic SOI, which couples spin mag- 
netic moments to the lattice, tends to restore the com- 
mensurability in the system. The magnetic moments lie 
mainly in the ab-plane (apart from small canting in the 
c-direction) [19j. The q y =l structure is mainly formed 
by the NN interactions and the magnetic anisotropy [12| , 
while the LR interactions are not affective due to the pe- 
riodicity constraint. For the q y = \ structure, the spins 
are aligned parallel to the longest Mn-0 bonds and min- 
imize the single- ion anisotropy energy. Moreover, the 
magnetic coupling along the b-axis minimizes the en- 
ergy of LR interactions. For the large- L structures, the 
situation is even more complex: certain spins minimize 
the single-ion anisotropy energy, while the directions of 
other spins compromise between NN, LR, and relativistic 
anisotropic and Dzyaloshinsky-Moriya interactions. For 
the q y = \ structure one can clearly distinguish different 



"domains" formed by the NN interactions and the single- 
ion anisotropy, respectively. Most importantly, all odd- 
periodic magnetic structures restore the inversion centers 
(associated with the central Mn-sites in Fig. [3]) and thus 
exclude any ferroelectric activity. 

Thus, we predict the oscillatory behavior of P depend- 
ing on q y =l/L (Fig. H]). The spin-spiral inhomogeneity 
in the even-periodic structures gives rise to the electronic 
polarization, similar to the E-state, which decreases with 
the increase of L. The canting of spins away from the 
collinear arrangement substantially reduces the values of 
P. For example, the electronic polarization |P|= 1.1 
/iC/cm 2 , obtained for the q y = z \ structure of TbMnOs 
(Fig. H]), is reduced by factor two in comparison with 
the collinear E-state (Fig. [2j). Taking into account the 
numerical values of <p (Fig. [3]), this reduction is readily 
explained by the angle-dependence of P in the nonrel- 
ativistic case (Fig. [2]). Thus, the SOI only forms the 
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inhomogeneous spin-spiral state, while the microscopic 
mechanism yielding finite P is essentially nonrelativistic. 

The existing theories overestimate the values of P 
(sometimes by an order of magnitude) @, [H| , which is 
one of the unresolved problems for the multiferroic man- 
ganies. In our calculations we obtain P||a= 1.4 /iC/cm 2 
for the q y = \ structure of HoMn0 3 (0=120°) and 0.3 
/iC/cm 2 for the q y = \ structure of TbMnOs shown in 
Fig. |4j Both values are about three times larger than the 
experimental ones [5[. Nevertheless, our theoretical ap- 
proach also suggests that it is difficult to obtain a good 
quantitative agreement because P is sensitive to many 
factors: (1) the precise value of (/>, which itself depends 
on the fragile balance of many magnetic interactions; (2) 
structural relaxation {zHHj; (3) due to the oscillatory be- 
havior (Fig. HJ, P can be reduced by possible deviations 
from the even-periodic commensurate alignment. 

Finally, we comment on the relative roles played by 
the magnetic and structural relaxation effects in the for- 
mation of inhomogeneous magnetic structures. The be- 
havior of the magnetic stabilization energies AE, defined 
as the energy difference between fully optimized magnetic 
structure with SOI and the homogeneous spin-spiral state 
without SOI, is explained in Fig. HI Since each increase 
of the supercell provides additional degrees of freedom 
for the magnetic relaxation, AE decreases with the in- 
crease of L. The absolute value \AE\~6 meV/Mn ob- 
tained for the q y =\ structure of HoMnOs is compara- 
ble with the energy gain caused by noncentrosymmetric 
atomic displacements in the E- state (~8 meV/Mn [181]). 
Thus, even for q y = \, the magnetic relaxation cannot be 
neglected and should be considered on an equal footing 
with the structural relaxation. Moreover, the transition 
to the inhomogeneous state is driven by the SOI, prior 
to structural relaxation, and in this sense there is no 
conceptual difference between HoMn03 and TbMnOs. 
Since \AE\ increases with L, the relative role of the mag- 
netic relaxation is also expected to increase. The mag- 
netic relaxation alone readily explains the experimentally 
observed q y ~\ periodicity in TbMnOs: AE of the rel- 
ativists origin shifts the total energy minimum of the 
homogeneous spin-spiral state (g y =0.68, Fig. [1]) towards 
smaller q y . Thus, the new theoretical minimum corre- 
sponds to q y =\ (Fig. HJ . Similar behavior was found for 
HoMnOs, which is expected to form similar fourfold pe- 
riodic structure. The disagreement with the experimen- 
tally observed twofold periodicity in HoMnOs is prob- 
ably caused by the neglect of the structural relaxation, 
which is relatively more important for q y —\- Quantita- 
tive aspects of interplay between magnetic and structural 
relaxation effects should be addressed in future theoreti- 
cal studies. At present, the structural relaxation cannot 
be easily implemented in the present model analysis. 

In summary, we propose that the homogeneous spin- 
spiral state in orthorhombic manganites is deformed by 
the relativistic SOI. The thus induced inhomogeneity is 



responsible for the ferroelectric activity of .RMnOs. It 
would be interesting to check our finding experimentally. 
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